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Abstract 



We perform a detailed study of the various one-loop contributions leading to produc- 
tion of the Standard Model Higgs boson in association with a hard jet. This production 
mode contributes to the current Tevatron exclusion limit of the Standard Model Higgs 
with 160 GeV < Mh < 170 GeV, and will also be important for discovery and interpre- 
tation of new scalar bosons at the Large Hadron Collider (LHC). We include top- and 
bottom-quark initiated contributions, maintaining the exact dependence on the quark 
masses, and also study previously neglected W- and Z-boson mediated effects which 
shift the qg and qq production modes. We consider the deviations from commonly used 
approximations for the Higgs boson transverse momentum spectrum caused by the fi- 
nite top-quark mass, bottom quark contributions, and electroweak gauge boson terms. 
All three effects act to decrease the Higgs boson transverse momentum distribution for 
observable momenta, with shifts reaching —8% at the Tevatron and —30% at the LHC. 
The shifts have a significant dependence on the Higgs pT, and are especially important 
if large momenta are selected by experimental cuts. 



1 Introduction 



The Higgs boson is the last undiscovered particle of the Standard Model (SM), and its 
discovery is a major goal of both the Tevatron and the Large Hadron Collider (LHC) physics 
programs. Direct searches at LEP restrict the SM Higgs boson to have a mass Mh > 114.4 
GeV [1]. Recently, the Tevatron collaborations have ruled out the mass range 160 GeV < 
Mh < 170 GeV at 95% C.L. [2]. If a scalar particle is discovered at either collider, the 
measurement of its properties will be crucial to determine whether the particle found is the 
Standard Model Higgs boson, or whether it hints at physics beyond the SM. 

The theoretical community has devoted significant effort to understanding precisely the 
production cross section and decay widths of the SM Higgs particle in order to facilitate 
such studies. The dominant production mode at both the Tevatron and the LHC is the 
partonic mechanism gg ^ H proceeding through a top-quark loop [3-7]. The exact next-to- 
leading order (NLO) QCD corrections to this process are known, including the corrections 
to the smaller bottom-quark induced contribution [8,9]. In the effective theory with the 
top quark integrated out by taking — > oo to produce an HGG operator, both the NLO 
corrections [10] and the NNLO corrections are known [11-13]. When normalized to the full 
mrdependent leading-order result, this effective theory reproduces the exact NLO result 
to better than 1% for < 2mt and to 10% or better for Higgs boson masses up to 1 
TeV [14]. The NNLO differential distributions in the effective theory were first obtained in 
Refs. [15-18], and detailed studies of the effects of experimental cuts on Higgs boson cross 
sections have been performed [19-22]. Resummation of logarithmically-enhanced threshold 
corrections to the cross section has been studied [23-25]. It was recently demonstrated that 
corrections arising from the Higgs coupling to W- and Z-bosons could affect the gg ^ H 
cross section at the 6-7% level [26-28]. The inclusion of such corrections is important for 
precision predictions of Higgs boson properties, and also impacts the exclusion limits set by 
the Tevatron collaborations. The theoretical uncertainty from missing higher-order terms is 
estimated to be roughly ±10% after these QCD and clcctroweak effects arc accounted for. 

Production of a Higgs boson in association with one or more jets is an important process 
for experimental searches. For example, the study of pp — > i?-|-jets with the subsequent decay 
H '-ff [29] or if — > WW [30] has significantly different kinematics than its background, 
and offers a good potential for Higgs discovery at the LHC. At the Tevatron, the multi- 
jet bins contribute roughly 50% of the exclusion limit [31]. Since the Higgs boson obtains 
a transverse momentum from recoil against hard jets, the study of Higgs+jet production 
is related to that of the Higgs transverse momentum spectrum. The calculation of Higgs 
boson production in association with a jet proceeding through top-quark loops, maintaining 
the exact dependence on the mass of the internal quark, was first presented in Rcf. [32] 
and further studied in Ref. [33]. The NLO QCD corrections to the Higgs boson transverse 
momentum spectrum have been computed in the effective theory with — > cxo [34-36], and 
the effect of bottom-quark induced contributions has been quantified [37] . The resummation 
of logarithmically enhanced corrections in the low pr limit have been considered [38-42]. 
In addition to aiding the discovery of the Higgs boson, study of its transverse momentum 
distribution has been discussed as a probe of physics beyond the SM [43-46] . 
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In light of the current experimental sensitivity at the Tevatron to Higgs production in 
association with a jet, and also in preparation for future studies at the LHC, we perform 
a detailed study of the various contributions to the one-loop Higgs transverse momentum 
spectrum in the SM. Some of these effects have previously been studied separately; we 
combine them into a single calculation and consider their collective impact. We study the 
effect of the top-quark mass and bottom quarks on the gg Hg, qg Hq, and qq Hg 
partonic processes. Since electroweak effects were shown to have a sizable effect on the 
inclusive cross section, we also consider W and Z boson mediated electroweak contributions 
to the Higgs spectrum. These affect both the qq and qg partonic processes, and were first 
considered in Ref . [47] . However, the interference of these such corrections with the top-quark 
initiated diagrams, which we find to be the dominant effect, was not studied. For observable 
values of the Higgs transverse momentum at both the Tevatron and the LHC, the effects of 
the finite top-quark mass, bottom quark terms, and electroweak terms all act destructively 
to reduce the high pj- cross section obtained from the — cx) effective theory. For relevant 
momenta at the Tevatron, the reduction reaches —8%. At the LHC, the reduction can reach 
—30% for Higgs transverse momenta with appreciable rates. We combine the exact t, 6, W , 
and Z mediated contributions in a FORTRAN code which allows arbitrary cuts on the Higgs 
boson kinematics to be imposed. 

Our paper is organized as follows. In Section [2] we present details of our analytic compu- 
tation of the various terms which give rise to the Higgs px distribution. We present numerical 
results for the Tevatron and for the LHC with y/s = 10 TeV in Section [3l We conclude in 
Section HI 

2 Calculational details 

We describe our calculation of the Higgs boson transverse momentum spectrum at one-loop 
accuracy, including terms induced by both heavy quarks and electroweak gauge bosons. We 
begin with the partonic channels affected by both quarks and electroweak corrections, the 
qq and qg initial states. Example diagrams contributing to the partonic processes qq Hg 
and qg — > Hq are shown in Fig. [H As the two processes are related by a crossing symmetry, 
we demonstrate the calculational details for the process q{pi) + q{p2) ^ H{j)h) + gips)- 
The amplitude for this process can be expanded in terms of two spin structures and their 
associated form factors, 

^ = T^ Yl ^''ii,u,s,MH,m,)ri + T^{u,i,s,MH,m,)r^„ (2.1) 

x=W,Z,b,t 

where we have defined the Mandelstam invariants t = {pi—ps)'^, u = {p2~P?)'^i s = (pi+P2)^- 
The sum over x includes all internal states with relevant couplings to the Higgs boson. We 
include W and Z bosons, top quarks, and bottom quarks. The spin structures can be written 
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as 

= Tlj- [pi ■ P3 v\il{v^ + 0^75)^1 - Pi ■ 63 v\i)'i{v^ + a^75)Mi] , 

= 7;;;. [p2 ■ P3 ^2/3(^.' + «x75)t^i - P2 ■ £3 ^2l^3(^x + ax.75)«i] • (2-2) 

We note that this process conserves chirahty, and is not suppressed by the initial-state hght- 
quark masses. The relevant couplings are simple to derive in terms of standard Feynman 
rules, and are found to be as follows: 

vw = gsgMw[iwiy + iwlY], aw = 2g,gMwwlwl, vz = gs—Mz [{zlf + [zlf] , 

az = 2g,^Mzzlzl vt,b = -glg^, a^,, = 0, = wt - ^ 



'cw " • ^^^2M^' ^'^ ' ^ ^ 2v^' 
'"'^ " """" 2V2' ^'cwU 3'^)' Acw' 



Here, gs and g are respectively the strong and weak coupling constants, and cw, sw are 
the cosine and sine of the weak mixing angle. We have parametrized the couplings of the 
quarks to the W, Z bosons as i {w^ + ^075) and i {z^ + ^075), respectively. For w^^ and z^ ^ 
the superscripts u and d denote respectively arbitrary up- and down-type quarks. 

One interesting fact to note about the amplitude defined in Eqs. fl2.ip and (12.21) is that 
it vanishes when becomes proportional to either pi or p2. As the electroweak corrections 
have no contribution from t-channel gluon exchange, this suggests that they will have a 
harder pt spectrum than the typical contributions, which is indeed what we find later. 
To derive expressions for the form factors JF^, we use standard techniques to reduce the 
contributing diagrams to a set of master integrals. We implement this reduction using the 
packages FORM [48] and QGRAF [49] . We give analytic expressions for the form factors in 
the Appendix. We note that the expressions for the top and bottom quark were first derived 
in Ref. [32], but we reproduce them here for completeness. 

H 





9 / \ 9 

Figure 1: Example one-loop diagrams contributing the the partonic process qq Hg 
through internal W, Z bosons (left figure) or fermions (right figure). 

Given the expression for the amplitude, it is straightforward to derive the differential 
cross section. In terms of the partonic center-of-momentum frame scattering angle 9 of the 
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final-state Higgs boson, the differential cross section is 

da 11(1- Mfj/s) 



d cos9 36 2s 



167r 



E 

spin,color 



M\ 



(2.4) 



with s + 1 + u = Mfj. The differential result for the qg channel can be obtained by crossing 
t ^ s and by changing the initial-state color averaging. 





H 



Figure 2: Example one-loop diagrams contributing the the partonic process gg 
through internal top and bottom quarks. 



Hg 



We also study the numerical effect of bottom quark and nit ^ oo corrections to the 
dominant partonic process leading to production of a Higgs boson plus one jet, g{pi) + 
9{P2) —>■ H + gips)- We have reproduced the cross section for this partonic process, which 
was first derived in Ref. [32]. Example diagrams are shown in Fig. [2l The amplitude for this 
process takes the form 



M 



99 



E 

x=b,t 



9s9 

IQ^2 ^ 2Mw 



Qi{i,u,s,MH,m^)Ti 



(2.5) 



with the same definition of Mandelstam invariants as utilized above for the other channels. 
We choose the following basis of spin structures to describe this process: 



Fi 



,a rabc 



When computing |A^p, we use the following sums over the various polarizations: 



J2e^,ie:r = -g'^^ + 



Px^x + Px^^ 
Px ' ^x 



pols 



(2.6) 



(2.7) 



with rii = P2, n2 = Pi, and = pi. The expressions for the form factors Qi are again 
straightforward to derive, and are given in the Appendix. The differential cross section for 
the gg partonic channel takes the form 



da 



1 1 (1 - M^/S) 



rfcos^ 256 2s 



167r 
4 



E 

spin,color 



99 \ 



(2.8) 



All amplitudes considered in this paper are finite and do not require renormalization. 
Wc have checked both analytically and numerically that our exact heavy-quark mediated 
contributions reproduce the effective theory results upon expansion in the fermion mass. We 
have written a FORTRAN code implementing the heavy-quark and electroweak contributions 
to Higgs+jet production valid for arbitrary Higgs boson mass, and which allows phase-space 
constraints to be imposed. We now use this code to study the numerical impact of the 
electroweak corrections, the bottom quark, and finite-mt effects, on the Higgs transverse 
momentum spectrum at the Tevatron and the LHC. 

3 Numerics 

We present numerical results to quantify the effect of various contributions to the Higgs 
transverse momentum spectrum. We utilize the parameter values that follow: 

Mw = 80.399 GeV, = 91.188 GeV, 
mt = 172.4 GeV, m5 = 4.8GeV, 

Gf = 1.16639 X 10"^ GeV-l (3.9) 

To define the remaining electroweak parameters we choose Mw, Mz, and as inputs 
and employ tree-level relations between the various quantities. We use MSTW 2008 parton 
distribution functions [50] extracted to NLO accuracy, and the consistent choice of with 
NLO running. For the studies we perform here we set = /ip = Mh- We note that 
the 6-quark contributions are small enough that our use of the pole mass versus a running 
mass has very httle effect on the numerical results. We set = 1.96 TeV for the Tevatron 
center-of-mass energy, and = 10 TeV for the LHC. To illustrate the sizes of the considered 
effects we set Mj^ = 120 GeV. 

We compute all cross sections at the one-loop level, which are the leading contributions to 
the Higgs pt spectrum. The next-to-leading order corrections to the top-quark contributions 
are known in the limit mt — > oo, but corrections to the W , Z, and bottom-quark terms we 
also consider are not known. Our primary results are the fractional deviations induced by 
electroweak effects, 6-quarks and the finite top-quark mass. If the corrections to the various 
terms are similar to that of the top quark, then the fractional deviations found here are 
unchanged. Since most of the effects are destructive interferences which reduce the rate, 
using the deviations found here leads to a conservative estimate if the corrections to the 
W , Z, and h terms are smaller. For the top quark, an alternative approach is to attempt to 
include the effective theory i^-factor to estimate the NLO corrections to thep^ spectrum [51]. 

3.1 Tevatron 

We begin by discussing the numerical effects of the various terms on the transverse momen- 
tum spectrum at the Tevatron. The dominant partonic process leading to production of a 
Higgs boson in association with a jet is gg — > Hg, which occurs through the Higgs coupling 
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to top and bottom quarks. The contributions from the qg channel contributes roughly half 
as much as gg while the qq channel is much smaller. For the gg and qg channels and when 
studying the fractional deviations from the effective theory result, we implement the cut 
Pt > 15 GeV, consistent with the CDF and DO definitions of a jet [2]. 



3.1.1 Effects of finite fermion masses 

A typical approximate method for computing this cross section involves integrating out the 
top quark in the limit 2mt ^ to produce the effective Hgg operator 



Ceft = -as^HG^.G'^^'' . (3.10) 



The cross section is computed using this Lagrangian and then normalized to the exact mi- 
dependent gg —y H interaction, leading to the approximate result 

(-.exact \ 
CT^^^. (3.11) 

We begin by quantifying how well this approximation works as a function of Higgs pt- We 
first show in Figs. [3] and H] the pt spectra for each partonic channel separately for three cases: 
the approximate cross section including just the top quark defined in Eq. (13. lip : the result 
with the exact top-quark mass dependence; and the cross section including both the top and 
bottom quarks with their exact mass dependence. These are shown to display the size of 
the cross section for various initial states and in each px bin. The sizes of the shifts are not 
easily seen on a logarithmic scale. We also define the deviation of the exact full cross section 
from the approximate form: 

^ 

^approx ' ^ ' 

where x indicates either the top-quark result or the cross section including both top and 
bottom quarks. Because the gluon parton distribution function rises rapidly at small Bjorken- 
x, most events in the qg and gg channels occur with small Mandelstam invariants. This 
implies that the loop integrals can be expanded around large m^, and the effective theory 
approximation works well. This is not true for the qq channel, where events with large s 
are not highly suppressed. The deviations from the approximate result are shown in the 
right panel of Fig. |H Bottom-quark loops in the gg and qg channels lead to a destructive 
interference effect of up to —8% at pt < 50 GeV. The effect of the bottom quarks has a 
dependence on the value of pt, varying from —8% at pr ~ 15 GeV to +4% at pr ~ 100 
GeV. At high pt above 100 GeV, where the relative size of the qq channel is increased, 
the large positive deviations caused by this channel reach +30% or more. However, the 
contribution of this kinematic region to the total rate is very small, as is clear from the 
plots. The kinematic feature near px ^ 150 GeV in the exact qq distribution occurs because 
the constraint > ^/ p^ + Mfj + px leads to > 2mt and the onset of an imaginary part 
in the amplitude at this value. 
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Figure 3: Transverse momentum spectra for the gg (left panel) and qg (right panel) channels 
at the Tevatron. A cut pt > 15 GeV is imposed, and a bin size of 10 GeV is used. 
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Figure 4: Transverse momentum spectrum for the qq partonic channel at the Tevatron (left 
panel). A bin size of 10 GeV is used. The right panel shows the deviations from the 
approximate cross section for the top quark with its exact mass dependence and for the top 
and bottom quarks. 
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3.1.2 Effects of electroweak corrections 



We now study the effects of including tfie W and Z boson corrections. Tliese affect only the 
qg and qq partonic channels. As noted before, these are expected to affect Higgs production 
predominantly at high pt- We show in Fig. [5] the pt spectra of both the qq and qg channels 
for the cases when the top and bottom quarks are included with their exact mass dependence, 
and when the W and Z are added. Although the effects on the qq channel are quite large, 
this channel is highly suppressed because of its lower luminosity. The addition of the W 
and Z-mediated terms contributes an additional destructive interference at high transverse 
momentum in the qg channel. 
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Figure 5: Transverse momentum spectra for the qq and qg channels at the Tevatron, with 
and without the electroweak corrections. The exact m/-dependent amplitudes have been 
used to obtain both distributions. 



The fractional deviations in each of the affected partonic channels and in the total cross 
section are shown in Fig. [6l The deviations are measured from the cross sections including 
both the top and bottom-quark terms with their exact mass dependence. Although the gg 
partonic channel is the largest by a factor of a few, the electroweak effects are large enough 
to result in a significant destructive interference for high Higgs px- The effects reach —4% 
for Pt ~ 50 GeV and —8% for px ~ 100 GeV. Also shown is the total deviation form the 
approximate top-quark result defined in Eq. 03.111) arising from W, Z, b loops and the exact 
top-quark mass dependence. The total deviation reaches —8% at pt ~ 15 GeV and remains 
at -3% from pr ^ 50 - 100 GeV. 



3.2 LHC 

We now repeat our numerical study for the LHC, assuming ^/s = 10 TeV. To avoid singu- 
larities and the effects of large logarithms, we implement the cut px > 30 GeV in the qg and 
gg partonic channels. 
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Electroweak Deviations Total Deviation 




Figure 6: Fractional deviations of the bin-integrated cross sections including the W and Z 
from that using the top and bottom quarks with their exact mass dependence (left panel) at 
the Tevatron. Results are shown for the qg and qq channels, and for the total cross section. 
Also shown is the total deviation from the approximate top quark result defined in Eq. (13.111) 
arising from W, Z, b loops and the exact top-quark mass dependence (right panel). 

3.2.1 Effects of finite fermion masses 

We again begin by studying the effect of using the approximate cross section defined in 
Eq. (13. lip versus the exact result using the full top and bottom quark mass dependence. 
The transverse momentum spectra for the various partonic channels are again shown in 
Figs. [7] and [8] for orientation. One important difference from the Tevatron case is that the 
large deviations in the qq channel do not become visible over the gg and qg channels even for 
transverse momenta of a few hundred GeV. The dominant effects up to pt = 300 GeV are 
therefore the deviations caused by using the exact mt-dependent amplitudes and including 
bottom quarks in the gluon-initiated channels. These are shown in Fig. [HI they can be 
quite large and have a complicated dependence on px due to cancellation between the shifts 
caused by the top and bottom-quark terms. The deviations resulting from including both 
the exact top and bottom-quark mass dependence are —4% for 30 GeV < Pt < 50 GeV, at 
the percent-level or less from 50 GeV until 150 GeV, and then become increasingly negative, 
reaching -30% at pr ~ 300 GeV. 

3.2.2 Effects of electroweak corrections 

We now study the effect of adding W and Z-induced electroweak corrections at the LHC. 
We expect them to be smaller than at the Tevatron, as the gluon-initiated partonic channels 
dominate to a larger extent at the LHC. The transverse momentum spectrum for the qq 
and qg channels, both with and without the electroweak terms, are shown in Fig. [91 Again, 
the deviations in the qq channel are large, but the numerical impact of this component is 
suppressed by its low luminosity. 
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Figure 7: Transverse momentum spectra for the gg (left panel) and qg (right panel) channels 
at the LHC. A cut pt > 30 GeV is imposed, and a bin size of 10 GeV is used. 




Figure 8: Transverse momentum spectrum for the qq partonic channel at the LHC (left 
panel). A bin size of 10 GeV is used. The right panel shows the deviations from the 
approximate cross section for the top quark with its exact mass dependence and for the top 
and bottom quarks. 
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Figure 9: Transverse momentum spectra for the qq and qg channels at the LHC, with and 
without the electroweak corrections. The exact mj-dependent amphtudes have been used to 
obtain both distributions. 



The fractional deviations for each partonic channel and for the total cross section caused 
by the electroweak effects are shown in Fig. (TDl As expected, the results are smaller than 
at the Tevatron. The W and Z-mediated contributions induce a —3% shift for pt values 
between 100 and 300 GeV. Also shown is the total deviation form the approximate top- 
quark result defined in Eq. ( 13. lip arising from W , Z, h loops and the exact top-quark mass 
dependence. This shift shrinks from —4% at px ~ 30 GeV to —1% at px ~ 100 GeV, and 
then continuously grows to reach —32% at pr = 300 GeV. 



4 Conclusions 

In this paper we have performed a detailed study of the possible one-loop contributions to 
the Higgs boson transverse momentum spectrum. We have considered effects arising from 
bottom quarks and the Z electroweak gauge bosons, as well as corrections arising from the 
full mass dependence of the top quark. The deviations from the often-used approximation 
defined in Eq. (13.111) . where only an infinitely heavy top quark is included, are generically 
negative, reaching —8% for moderate pt at the Tevatron and —30% for transverse momenta 
of several hundred GeV at the LHC. The corrections have a non-neglible dependence on the 
Higgs Pt- At the LHC, they may be particularly important if regions of large transverse 
momenta are selected, which has been discussed to increase to discovery potential in several 
Higgs boson search channels [29,30]. The electroweak corrections reach —8% at the Tevatron, 
and could potentially have an effect on the current exclusion limits. As far as we are aware, 
these contributions have not been implemented in any available calculation, nor in any 
experimental analysis. If a scalar is observed at the LHC, these effects should be accounted 
for to avoid mistakenly claiming a hint of new physics. 



11 




Figure 10: Fractional deviations of the bin-integrated cross sections including the W and Z 
from that using the top and bottom quarks with their exact mass dependence (left panel) at 
the LHC Results are shown for the qg and qq channels, and for the total cross section. Also 
shown is the total deviation from the approximate top quark result defined in Eq. (13.111) 
arising from W, Z, b loops and the exact top-quark mass dependence (right panel). 
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Appendix 

We present here the form factors that appear in the various partonic processes. We begin 
with the form factors jF^(t, m, s, Mj;^, m^) defined in Eq. (12. ip that compose the qq Hg 
and qg —>■ Hq channels. We recall that the Mandelstam invariants satisfy s + t + u = Mfj. 
For the fermionic contributions, where x = t,b, we find 



Sruf rBubl(M^,m^) - Bubl(s,m}) 



s-Mjj\ s-Ml 




(5.13) 
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For the electroweak contributions with x — W,Z,we find 

(-£ + M^) {-iu + M^u + i) Box(t, M, , M^) 



i (i+s) 



^ Bub2 (^,M2) ^ + (-£+ M^) Tri (^^ M^, M^) Bubl(M^,Mi 
"2.^~l~2 . ~l~4 

2 (-s^f + s^M^ + 2M^si- £2 ^ £3^ ^2 ^ ^2 ) 

+ s) s P 

J-£+M^)^Bub2(n,M2) H+M^)ln(l-^) 



-?2 + M^) s ^2 



5^2 



+2 {-3Pu-5ius + 2M^si + 3suM^ + 3M^iu)\n (l - (f+s) ^ 



xs^^r^u"^ + 8 



{-u + M^)sP {i+s) ' 



(5.14) 



The various auxihary functions that appear in this expression are the finite parts of loop 
integrals that appear in the calculation; all poles that appear in each integral cancel in the 
amplitude. They are given by finite parts of the following expressions: 



Bubl(p^m^) = 

Bub2(p^m2) = 
TVi(p2,p2 r„2) ^ 

TViF(p2,p2 r„2) ^ 

Box(p^,p2^p2^77i2^ = 
X 



1 



^2 _ j^2 _|_ pj2 

d'^k 1 1 



l'jld/2 f^2 (J^ _|_ _ ^2 

d'^k 1 1 



finite 



finite 



1 



with {pi 

d'^k 



J^2 (J^ _|_ _ ^2 (J^ p^j^ p^-j2 

' = 0, 



finite 



P2) 



1 



1 



j^d/2 1^2 _ ^2 (^/j, _|_ p^'j2 _ ^2 (^/j. _|_ _|_ p2)2 — 777,2 

with (pi +^2)^ = 0, 

d'^k 1 1 1 

77r<i/2 ^2 _|_ p^y _ ^2 _|_ p2 — ^3)2 _ 7772 
1 



{k+pi +P2-P3y 



with (pi - Pa) = 0, 



finite 



(P2 - Psf = 0, (pi + P2 - pa)^ = 0. 



(5.15) 



It is straightforward to express these integrals in terms of logarithms and dilogarithms using 
standard techniques [52]. All integrals were checked numerically against FF [53] for a variety 
of phase-space points. A further useful check on the exact calculation is comparison with 
the approximate forms of the amplitude in the limit that the internal loop particles become 
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much heavier than the Higgs boson or the Mandelstam invariants. By direct expansion of 
the relevant diagrams, we can derive the following simple expansions for the form factors: 

... A 



Srrif s 

T''ii,u,s,MH,Mv) """^"^ (5.16) 

We note that these expanded forms do not provide a suitable numerical approximation of 
the electroweak effects. 

For the Qi form factors defined in Eq. (12.51) . we first note that only two are independent 
due to the following symmetry relations: 

gi(t,u,s,MH,mf) = g3{s,u,i,MH,mf), 

gi{i,u,s,MH,mf) = -g2{u,s,i,MH,mf), (5.17) 

which imply that only and ^4 are needed to describe the amplitude. These take the 
following forms: 

mr i (Snif^u - sii - + 2 i s) BoxT (s, i, M|, m/^) 
gi{t,u,s,MH,mf) = ^ 

ruf {-i s + 8 m/^t - t'^ + 2 s u) BoxT (5, u, Mj, m/^) ^ ruf i (s - u) Bubl (t, m/^) 



i [u + s) s u 



nif i (4 m/s + i u) BoxT [i, u, M|, m/^) (^2 + p + 4 t Bubl (s, mj 



2\ 



+2 — ^ - ■ - ^ — - 4 

^ rrif {s - i) Bubl (n, mj^) ^ m/ PTriF2 (f, m/^) (s^ - ^2^) ^ m/ u TriF2 (n, m/^) (s^ - t^) 

s [i + s) u^s"^ Pi 

+4 [8 Pu + Pus + 2 + i^s - i^u + 8 s^i + 4iu - u^i - 3 u^P - 3 u^P 

+2 fu^ + s u^i + s^u^ + s^P + su'^ + 2s u'^P] mj Bubl (M|, mj^) 

r , , .,2, ^1^1 rrif 5 TriF2 (s, mf^) (? + u^) (su + is + iu)'mf 

x\su (t + s) (u + t) (u + s)\ -2 ^ 1 ^ ^ ^ ^ - 8 ^ 7 ^ 

\ ^ ^ ^ ' ^ tu^ {u + t)su 

+2 ruf {Pu s + Ps + 4 i mj'^iP - 2 u^P - 2 tPP + s u^i + 16 PiPnif'^ 

+4 Pm/^u + 5 n^) TriF (s, M|, m/^) { t (u + Pj sipy^ + 2 {s^i - s^u + iu 

-tPs^ + Anif'^u f - tips - Am/iPs - t vP) nif TriF (f, M|, my^) u^^s"^ 

(^^^^^ + s^?2 + t + 4 nif'^Pi - PP - Pus - A rUf'^Ps - Pu) nij TriF {u, Mj, m/) _ 
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ruf s (12 m/u - + 4 f s) BoxT (g, f, M|, m/^) 

g4{t,u,s,MH,nif) = -2 

m/ s (4 m/^ - u) BoxT (s, u, M|, m/^) m/ (4 m/^ - u) BoxT u, M|, m/^) 

t u u 
{t + 2u) nif sBubl (s, m./^) s {2u + s) rrif Bubl (t, m/^) nif {t s — u^) 



-16 ^ ' — ^ ' ^ ' _ ^ / ' ^ „ ^ ' ' ' - 16 



+16 

+4 



(t2 + U2 (U + 5)^ U2 ui{u + S) [u + 

m/ s (4 ^ i2 s + 2 P?2 + s'^i + s + 2 s'^u + b s v? + bi + A u^) Bubl (M|, m/^) 
m/ (4 s + 2 i^s + 2 f u^s + u"^ + 4 m/2f2^i _ f2^2 _ 4 ^^2^3 _^ g ^^2^2^^ 



2-3 



xTriF (s, M|, m/^) + 4 m/ (2 s^i + 4 f u + 2 £ u^s - ?2^s^ + + 4 m/^u - 4 mfv: 
+8 TriF (i, M|, m/) { (n + ,s) t {L^y' - 4 ^^'^ ""^^ " 



t u 



2 



xMF(«.Mg. V) -S"^-"^"^?''-"^') -8'"^'^'™^('''"^T (5.19) 



Two new integrals appear in these form factors: 
TWF2(p;,m=) = / 



BoxT(pi,p2,P3>"^ ) 



rn? (k + pi)2 — 777,2 _|_ _|_ ^2)2 — 777,2 
with (pi+p2)' = 0, p^ = 0, 

rf^A; 1 1 1 

i7r'^/2 /j2 _ ^2 ^j^ _|_ p2^2 _ ^2 _|_ — p3)2 — 7772 

X 77— ^ r- ^ with {pi - = 0, 

{k + pi+p2- P3> - 

(P2-P3)' = 0, (Pl+P2-P3)'=0. (5.20) 

It is again useful to study the approximate expressions for these form factors as 777/ becomes 
large. They take the form 

gi{t,u,s,MH,mf) — , 

Srufsu 

g4{i,u,s,MH,mf) ""'^^ (5.21) 

37?7/t 
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